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Abstract — Multi-cell cooperation (MCC) is an approach for 
mitigating inter-cell interference in dense cellular networks. 
Existing studies on MCC performance typically rely on either 
over-simplified Wyner-type models or complex system-level sim- 
ulations. The promising theoretical results (typically using Wyner 
models) seem to materialize neither in complex simulations nor 
in practice. To more accurately investigate the theoretical perfor- 
mance of MCC, this paper models an entire plane of interfering 
cells as a Poisson random tessellation. The base stations (BSs) 
are then clustered using a regular lattice, whereby BSs in the 
same cluster mitigate mutual interference by beamforming with 
perfect channel state information. Techniques from stochastic 
geometry and large-deviation theory are applied to analyze the 
outage probability as a function of the mobile locations, scattering 
environment and the average number of cooperating BSs per 
cluster, L For mobiles near the centers of BS clusters, it is shown 
that outage probability diminishes as 0(e~ e 1 ) with < v\ < 1 
if scattering is sparse, and as 0(£ ~" 2 ) with V2 proportional to 
the signal diversity order if scattering is rich. For randomly 
located mobiles, regardless of scattering, outage probability is 
shown to scale as 0(£" 3 ) with < v-j, < 0.5. These results 
confirm analytically that cluster-edge mobiles are the bottleneck 
for network coverage and provide a plausible analytic framework 
for more realistic analysis of other multi-cell techniques. 



I. Introduction 

Inter-cell interference limits the performance of cellular 
downlink networks but can be suppressed by multi-cell coop- 
eration (MCC). The existing high-speed backhaul links allow 
base stations (BSs) to exchange data and channel state infor- 
mation (CSI). Thereby, cells can be grouped into finite clusters 
and BSs in a same cluster cooperate to decouple the assigned 
mobiles [l]-[4]. Despite extensive research conducted on 
MCC, the fundamental limits of cellular networks with MCC 
remain largely unknown due to the lack of an accurate and 
yet tractable network model. This paper addresses this issue 
by proposing a novel model constructed using a Poisson point 
process (PPP) for BSs and a hexagonal lattice for clustering 
said BSs. Based on this model, techniques from stochastic 
geometry and large-deviation theory are applied to quantify 
the relation between network coverage and the average number 
of cooperating BSs. 
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A. Modeling Multi-Cell Cooperation 

Quantifying the performance gain by MCC requires accu- 
rately modeling the cellular-network architecture and account- 
ing for the relative locations of BSs and mobiles. These factors 
are barely modeled in Wyner-type models where base stations 
are arranged in a line or circle, interference exists only between 
neighboring cells and path loss is represented by a fixed 
scaling factor [5]. Due to their tractability, Wyner-type models 
are commonly used in information-theoretic studies of MCC 
[4], [6], [7], but fail to account for mobiles' random locations 
[8] and finite BS clusters in practice due to a constraint on the 
cooperation overhead [1], [9], [10]. The traditional hexagonal- 
grid model provides a better approximation of a practical 
cellular network, however, at the cost of tractability [11]. An 
alternative modeling approach is to model BSs using a PPP 
and construct cells as a random spatial tessellation [12]. The 
random model captures cell irregularity, is about as accurate as 
the hexagonal-grid model, and allows analysis using stochastic 
geometry [13], [14]. 

Building on [12] which assumes single-cell transmission, in 
this paper BSs are modeled as a homogeneous PPP that parti- 
tions the horizontal plane into Voronoi cells. Mobiles in each 
cell are randomly located and time share the corresponding 
BS. BSs are then clustered using a larger hexagonal lattice 1 to 
cooperate by interference coordination where BSs in the same 
cluster mitigate interference to each others' mobiles by zero- 
forcing beamforming that also achieves transmit-diversity gain 
[15]. Furthermore, to cope with fading, channel inversion is 
applied such that received signal power is fixed. This scheme 
is considered for simplifying analysis and can be implemented 
in practice by combining a transmit-diversity technique and 
automatic gain control widely used in code-division-multiple- 
access systems. It is worth mentioning that channel inversion is 
found in this research to reduce outage probability compared 
with fixed-power transmission. Outage probability specifies 
the fraction of mobiles outside network coverage for a target 
signal-to-interference ratio (SIR), assuming an interference 
limited network. This is the case of interest for MCC and of 
operational relevance for cellular networks. Let the average 
number of BSs in a cluster be denoted as £, called the 
expected BS-cluster size. This paper focuses on quantifying 

1 The hexagonal lattice is chosen arbitrarily for exposition. It is straight- 
forward to extend the current analysis to BS clustering using other types of 
regular lattice or random spatial tessellations by modifying the definitions of 
the variables p, p and D (defined in the sequel) based on the cell geometry. 
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the asymptotic rate at which outage probability diminishes as 
£ increases. 

This and any other clustering methods with finite cluster 
sizes and only intra-cluster cooperation have the drawback 
of cluster-edge mobiles exposed to strong inter-cluster inter- 
ference as quantified in the subsequent analysis. Intuitively, 
a better approach is to allow overlapping BS clusters for 
protecting cluster-edge mobiles. BS cooperation based on 
this approach can be implemented efficiently using belief 
propagation and message passing [16]— [18] but will eventually 
involve all BSs in the network and cause potential issues in- 
cluding overwhelming backhaul overhead, excessive delay and 
network instability. For these reasons, BS clusters in practice 
are usually disjoint [19]. This investigation suggests a much 
simpler approach for suppressing inter-cluster interference for 
cluster-edge mobiles by combining the current method of BS 
clustering with fractional frequency reuse [20] along cluster 
edges as discussed in the sequel. 

There exists a rich literature on analyzing outage probability 
for wireless networks with Poisson distributed transmitters 
[21]-[24]. Given that outage probability has no closed-form 
expressions [25], [26], a common analytical approach is to de- 
rive bounds on outage probability using probabilistic inequali- 
ties [27], which are sufficiently simple and tight for evaluating 
network performance given specific transmission techniques 
e.g., bandwidth partitioning [28] and multi-antenna techniques 
[29], [30]. The accuracy of these outage-probability bounds 
requires the presence of strong interferers for mobiles. Similar 
bounds for cellular networks with MCC can be loose since 
interference is suppressed using MCC. Therefore, this work 
deploys an alternative approach where large-deviation theory 
[31] is applied to quantify the exponential decay of outage 
probability as £ -> oo. A similar approach was applied in [32] 
to analyze the tail probability of interference in a wireless ad 
hoc network. 

B. Summary of Contributions and Organization 

To apply techniques from large-deviation theory, a new 
performance metric called the outage-probability exponent 
(OPE) is defined as follows. Since the network is interference 
limited and hence noise is negligible, the outage probability 
for an arbitrary mobile, denoted a.s Pq\jx, is given a.s 

Pout = Pr (y < 0) (1) 

where lo and I represent the fixed received signal power and 
random interference power, respectively, and 6 > is the 
outage threshold. Then the OPE is defined as 

= -log Pout ( 2 ) 
= -logPr(l>|) (3) 

where P out and I are functions of £ with £ omitted for ease of 
notation. It follows that deriving the scaling of ip{£) as £ — > oo 
yields the exponential decay rate of P out . Using large-deviation 
theory, simple OPE scalings are derived for different network 
configurations based on the rates at which the tail probabilities 
of random network parameters diminish as I — > oo. 



With interference being suppressed by increasing £, the 
network will eventually operate in the noise limited regime, for 
which the outage-probability for a typical mobile is either zero 
or one depending on if the received signal-to-noise ratio uj/cr 2 
is below or above 6. The value of lj depends on the average 
transmission power of BSs and channel distribution [see (10) 
in the sequel]. Therefore, the OPE becomes irrelevant for the 
case of a noise-limited network with channel inversion. 

The main contributions of this paper are summarized as 
follows. 

1) Consider a mobile located at the center of an arbitrary 
BS cluster, called a cluster-center mobile, and sparse 
scattering where beams have bounded amplitudes. Given 
MCC, the OPE for a cluster-center mobile, denoted as 
ip cc , is shown to scale 2 as follows: 

a) for the path-loss exponent a > 4, 

ci£ r< <p cc (£) r< -^t, £^oo; 

b) for 2 < a < 4, 

where C\ and c 2 are constants. 
This result shows that outage probability diminishes 
exponentially as £ — > oo for a high level of spatial 
separation (a > 4) or at least sub-exponentially if the 
level is moderate-to-low (2 < a < 4). 

2) Consider a mobile with a randomly distributed location, 
called a typical mobile, 3 and also MCC with sparse 
scattering. The scaling of the corresponding OPE is 
proved to be 

^(l-^jlog£^<fi(£)^^log£, £^oo. (4) 

This result implies that outage probability decays as 
I — » oo following a power law with an exponent smaller 
than 0.5. This decay rate is much slower than the sub- 
exponential (up to exponential) rate for a cluster-center 
mobile. The reason is that a typical mobile may lie near a 
cluster edge and consequently is exposed to strong inter- 
cluster interference. Comparing the outage-probability 
decay rates for cluster-center and typical mobiles sug- 
gests that cluster-edge mobiles are the bottleneck of 
network coverage even with MCC and protecting them 
from inter-cluster interference (e.g., assigning dedicated 
frequency channels) can significantly improve network 
coverage. 

3) Consider MCC with rich scattering modeled as Rayleigh 
fading. Note that fading affects the interference distri- 
bution but not received signal power that is fixed given 

2 Two functions f(z) and g(z) are asymptotically equivalent if — ¥ 1 

as ^ oo, denoted as /(*) ~ g(z); the cases of lim,^ > 1 and 

lim^oo < 1 are represented by f(z) > g(z) and f(z) -< g(z), 

respectively. 

3 A typical point of a random point process is chosen from the process by 
uniform sampling such that all points are selected with equal probability. 
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TABLE I 
Summary of notation 



Symbol 


Meaning 


(n , n cc 

r ' T 


OPE tor a (typical, cluster-center) mobile 


I, I cc 


Received interference power for a (typical, cluster-center) mobile 


I 


Expected BS-cluster size 


M 


Number of BSs in a typical cluster 


A 


PPP of BSs, density of $ 


n 


Hexagonal lattice for clustering BSs 


Y*,T*,U* 


Typical BS, BS-cluster center and mobile 


w 


Cluster of mobiles served by the typical BS cluster 


C(T r) 


Hexagon centered at T £ and having the distance v from T to the boundary 


/? /> 

P' p 


Distance from the center of a cluster region to an (edge, vertex) 




Mobile served by BS Y 


Ly 


Distance from BS Y to the affiliated mobile 




Transmission power for BS Y 


fy 


Beamformer used at BS Y 


h[/y 


Vector channel from BS Y to mobile U 


Q 


Path-loss exponent 




Outage threshold 




Fixed received signal power at a mobile 


AT, 


Signal diversity order for a typical mobile, the minimum value of N 


D 


Distance from a typical mobile to the boundary of the corresponding cluster 



channel inversion. The OPE for a cluster-center mobile 
is shown to satisfy 

(^av-l^\ogl^ip cc {C) r< \av\ogl, i -> oo 

where v > 1 is the minimum signal diversity order over 
different cells. It follows that outage probability decays 
as I — > oo following a power law with an exponent 
approximately proportional to a and v. By comparing 
the outage-probability decay rates for sparse and rich 
scattering, it is found that additional randomness in 
interference due to fading degrades the reliability of 
communications near cluster centers significantly. 
4) Last, the OPE scaling for a typical mobile with sparse 
scattering from (4) is shown to also hold for a typical 
mobile with rich scattering. The OPE scaling is largely 
determined by the probability that the mobile lies near 
cluster boundaries and outside network coverage due to 
strong inter-cluster interference. As a result, the scaling 
is insensitive to if fading is present, which, however, 
impacts the OPE scaling for a cluster-center mobile. 

The remainder of the paper is organized as follows. The 
network model is described in Section II. The OPEs with 
sparse scattering and with rich scattering are analyzed in 
Section III and Section IV, respectively. Simulation results 
are presented in Section V followed by concluding remarks in 
Section VI. The appendix contains the proofs of lemmas. 

C. Notation 

The complement of a set X is represented by X. The 
operator \X\ on X gives its cardinality if A" is a set or the 
distance from X to the origin if X represents a point in 
the plane R 2 . The superscripts T and f represent the matrix 
transpose and Hermitian transpose operations, respectively. 

The families of distributions having regularly varying and 
Weibull-like tails are represented respectively by RV(t) and 



WE(t) where r > is the index, and defined as follows. 
Define the distribution functions T and T of a random variable 
(rv) X as T{x) =_ Pr(X < x) and T = Pr(A > x). The rv 
X e RV(t) if F{x) = x~ T V(x) as x ->• oo with V(x) 
being a slowly varying function, namely lim ;E _ i . 00 = 1 

for all t > [33]. If X e WE(r), X has support [0, oo) and 
F(x) ~ cx^e-^ form some /3 e R, r e (0, 1), C > and a 
constant c> [34]. 

Other notation is summarized in Table I. 

II. Network Model 

A. Network Architecture 

The BSs are modeled as a homogeneous PPP $ = {Y} 
in the horizontal plane R 2 with density A where Y € R 2 is 
the coordinates of the corresponding BS. The mobiles form a 
homogeneous point process independent with $. By assigning 
mobiles to their nearest BSs, the horizontal plane is partitioned 
into Voronoi cells as illustrated in Fig. 1. It is assumed that the 
mobile density is much larger than the BS density such that 
each cell contains at least one mobile almost surely. Each BS 
Y serves a single mobile at a time, denoted as u(Y), selected 
from mobiles in the corresponding cell by uniform sampling. 
Consequently, the distance between an arbitrary BS Y e $ 
to the intended mobile, denoted as Ly, has the following 
distribution function [12]: 

Pr(Ly- > x) = e-* Xx \ x>0. (5) 

BSs are clustered using a hexagonal lattice VL = {T} where 
T e R 2 denotes the coordinates of a lattice point. Using 
the lattice points as cluster centers, the horizontal plane is 
partitioned into hexagonal cluster regions as illustrated in 
Fig. 1. Let C(T, r) denote a hexagon centered at T <G R 2 and 
having the distance r from T to the boundary. Thus the cluster 
region centered atTeO can be represented by C(T, p) where 
p is specified in Fig. 1. Note that p determines the density 
of the lattice The area of C(T, p) is 2\* r 3p 2 and hence 
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the expected BS-cluster size is I = 2^/3p 2 X. Let Y* denote 
a typical point in $, called the typical BS, and the mobile 
served by Y* is called the typical mobile and represented by 
U*. Moreover, define the typical cluster center T* E fl as 
one such that C(T*,p) contains Y* . The cluster of BSs lying 
in C(T*,p), namely $nC(r,p), is called the fy/"'cfl/ BS 
cluster, the associated cluster of mobiles is represented by 
W = {u(Y) | Y E $nC(T,p)}. 

B. Multi-Cell Transmission 

The cooperation in a BS cluster is realized using a practi- 
cal interference-coordination approach that requires no inter- 
cell data exchange [15]. Consider the typical BS cluster 
$DC(T*,p) and the affiliated cluster of mobiles U* . Assume 
that each BS employs Q antennas and mobiles have single- 
antennas. Let M denote the number of BSs and hence M = 
|$ n C(T*,p)\ is a Poisson random variable (rv) with mean 
£. It is assumed that Q > M so that each BS has sufficient 
antennas for suppressing interference to mobiles served by 
other cooperating BSs. As a result, Q is a rv and varies over 
different clusters. The analysis in the sequel focuses on the 
regime of a large average cluster size (£ — > oo) corresponding 
to the regime of large-scale antenna arrays (E[Q] — > oo). 
With expected deployment of large-scale arrays in future 
wireless networks [35], such an assumption may be viable. 
Furthermore, the analytical results will be shown to also be 
accurate for moderate numbers of antennas. For instance, it 
will be observed subsequently from simulation results (see 
Fig. 4) that for sparse scattering the derived asymptotic bounds 
on the OPE are tight for t smaller than 6 and E[Q] equal 
to I plus several more antennas to achieve moderate array 
gain. Let h$ Y E C represent the coefficient of the scalar 
channel from the fc-th antenna at Y to U and define the 
channel vector hj/y = [h\j Y , hy^ Y , • • • , h!ff Y ] for given Q. 
Moreover, let { Y E <C Q with ||fy|| = 1 denote the unitary 
transmit beamformer used at Y. The interference avoidance at 
Y* is achieved by choosing f Y to be orthogonal to the (M— 1) 
interference channels and the remaining N = Q — M + 1 
degrees of freedom (DoF), called the diversity order, are 
applied to attain diversity gain [36]. It is assumed that N > v 
with v > 1 being the minimum diversity order over different 
cells, where the constraint v > 1 ensures finite average 
transmission power under channel inversion for the case of rich 
scattering. Assuming perfect CSI at BSs, their beamformers 
are designed using the zero-forcing criterion as follows. 

Definition 1 (Interference coordination). Conditioned on Q = 
m, the beamformer fy» used at the typical BS Y* solves: 

maximize: |Fh[/»y»| 

subject to: tfh UY , = V U E U*\{U*} (6) 
f E C"\||f|| = 1. 

This algorithm is also considered in [37] for mitigating 
inter-cell interference in a two-cell network. Note that the 
computation of fy* requires Y* to have CSI of both the data 
channel and the (M — 1) channels from Y* to mobiles served 
by other cooperating BSs, which can be acquired by CSI 



feedback [38]. Given that the network is interference limited, 
with the beamformer designed as in Definition 1, the signal y 
received at U* is given as 



/ P Y 'f Y ,h U , Y ,X U ,+ Y \fPY^ Y ^U'YX u ( Y ) (7) 



Ye$nC(T*,p) 



where P Y denotes the transmission power of BS Y and xjj is 
a data symbol with unit variance and intended for U. Let S 
and I represent the signal and interference powers measured 
at U*, respectively. It follows from (7) that 



S = Py 



\f Y ,h. 



U*Y* 



and 



YE$nC(T*,p) 



P Y \f Y hu, Y \ 



(8) 

Besides mitigating interference using MCC, channel inver- 
sion is applied at BSs to cope with data-link fading. The 
transmission power P Y of BS Y is chosen such that the signal 
power received by the intended mobile is a constant w > 0. 
Consequently, S = lj and 



P Y » = 



LJ 



(9) 



where lj satisfies the average power constraint E[Py»] < P 
with P > and hence is given as 



P 



(10) 



U*Y* 



It is found in this research that channel inversion increases 
OPE (reduces outage probability) compared with fixed-power 
transmission. The reason is that fixed-power transmission 
causes fluctuation in received signal power that increases 
outage probability but can be removed by channel inversion. 
The analysis for the scenario of fixed-power transmission is 
omitted to keep the exposition precise. 

C. Channel Models 

The scattering environment affects the interference distribu- 
tion and hence the OPE. For this reason, both sparse and rich 
scattering are considered in the OPE analysis and their models 
are described as follows. 

1) Sparse Scattering: In an environment with sparse scat- 
terers, there usually exists a line-of-sight path between a 
transmitter and a receiver and fading is negligible compared 
with this direct path. Using beamforming in Definition 1, each 
multi-antenna BS forms a physical beam such that the main 
lobe is steered towards the intended mobile, nulls towards 
mobiles served by cooperating BSs, and side-lobes towards 
others [39]. This can be modeled such that the interference 
power I in (8) and transmission power Py» in (9) are given 
as 



^2 PyGu, y \Y 

YE$nC(T*,p) 



u* 



P Y , =ljLZ,WZ 



(11) 
(12) 



where the path-loss exponent a > 2, W Y is the main-lobe 
response of beamforming at Y, and Gjj Y is its side-lobe 
response in the direction from Y to U. In practice, the values 
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Fig. 1. (Left) The topology of the cellular network with Poisson distributed BSs clustered using a hexagonal lattice. The cells are drawn using thin lines 
and the cluster regions thick lines; BSs are marked using black dots. (Right) A hexagonal cluster region where p and p denote the distances from the cluster 
center to an edge and a vertex, respectively, and p = ^rp~. The cluster area is 2\/2tp 2 and hence the expected BS-cluster size is t = 2 v / 3p 2 A. 



of Wy and Gjj*y depend on the size and configuration of 
BS antenna arrays as well as transmission directions [39]. 
They are modeled as random variables (rvs) with the following 
properties. 

Assumption 1 (Sparse-scattering model). The rv Wy has 
bounded support [6, 6'] with 6' > 6 > 0. For U and 
Y associated with different BS clusters, the rv Gjjy has 
bounded compact support [0, 7] with 7 > 0. The set of rvs 
{Guy I U G U*,Y e $ nC(T*,p)} are independent and 
identically distributed (i.i.d.). 

For clarification, the equality 7 = 6' holds in theory since it 
is possible for a transmitter to direct a beam towards both an 
intended and an unintended receivers if they lie in the same 
direction. Nevertheless, given sufficiently sharp beams and 
randomly located nodes, such an event occurs with negligible 
probability and hence it can be assumed that 7 <C 6,6'. This 
assumption, however, is not required for the current analysis. 

2) Rich Scattering: The channel is assumed to be frequency 
non-selective and follows independent block fading. Rich 
scattering is modeled by i.i.d. Rayleigh fading as follows. 

Assumption 2 (Rich-scattering model). An arbitrary channel 
coefficient hfj Y is given as bJ£ Y = B^ Y \U — Y\~ a where 
B$ Y is a CAf(0, 1) rv. Any two rvs B$ Y and B [ ^', ] Y , with 
(k, U, Y) ^ (k', U', Y') are independent. 

It follows from Assumption 2 that an arbitrary channel 
vector hjjY can be written as hr/y = qjjy\U — Y\~ a where 
qjjy is a Q x 1 random vector comprising i.i.d. CJ\f(0, 1) 
elements. Moreover, the sequence {qr/y} is i.i.d. The signal 
and interference powers measured at U* are given by (11) 
and (12) but with the parameters Wy and Gjj*x re-defined 
as Wy* = \f Y *(lu*Y*\ 2 and Gu*y = |fyqc/*y| 2 - The lemma 
below follows from [36, Lemma 1] that studies zero-forcing 



beamforming (see Definition 1) for mobile ad hoc networks. 

Lemma 1 ([36]). For rich scattering and conditioned on N = 
n, Wy* is a chi-square rv with 2n DoF and {Gjj*y | Y € 
<J> <1C(T* , p)} are i.i.d. exponential rvs with unit mean. 

III. OPE with Sparse Scattering 

In this section, the OPE is analyzed for the environment 
of sparse scattering. Specifically, the OPE is characterized for 
a cluster-center mobile and for a typical mobile separately. 
The results show that mobiles near cluster edges limit network 
coverage. 

A. OPE for Cluster-Center Mobiles 

Consider a mobile located at the typical cluster center 
T* that is farthest from the interference zone among all 
mobiles and hence has the smallest outage probability, where 
an interference zone for a mobile refers to a region in the 
horizontal plane comprising interfering BSs. The OPE for 
a cluster-center mobile, denoted as (p cc , can be written by 
modifying (3) to account for the constraint U* = T*: 

ip cc (l) = -logPr(/ cc > e^u \U* =T*) (13) 

where I cc represents the interference power measured at T*. 
Asymptotic bounds on ip cc for large I are derived in the sub- 
sections and then combined to give the main result of this 
section. 

1) Asymptotic Lower Bound on the OPE: First, a lower 
bound on ip cc is obtained as follows. Slightly abusing notation, 
let T* also represent the typical cluster-center mobile. As 
illustrated in Fig. 2(a), C(T*,p) is the complete interference 
zone for T* . Therefore, I cc can be obtained by modifying (11) 
as 

I cc = P Y G t *y\Y -T*\- a (14) 

Ye*nC(T*,p) 
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(a) Cluster-center mobile (b) Typical mobile 



Fig. 2. The geometry of the interference zones (shaded regions) for (a) a cluster-center mobile and (b) a typical mobile, where the interferers are Poisson 
distributed with density A. The symbols T*, C(T*,p), Y*, and U* represent the typical cluster center, cluster region, BS and mobile, respectively, and 
L=\Y*-U*\. 



which is a power-law-shot-noise process [26]. It can be 
observed from (13) that the OPE is determined by the tail 
probability of I cc that, however, has no closed-form expression 
[26]. For the current analysis, it suffices by deriving an upper 
bound on I cc . This relies on decomposing I cc into a series 
of compound Poisson rvs inspired by the approach in [32]. 
To this end, the interference zone C(T*,p) is partitioned 
into a sequence of disjoint hexagonal rings {A n }^Li with 
A n = C{T* ,y/n + lp)\C(T* ,y/np). Note that {A n } have 
the same area as C(T*,p). The interference power measured 
at T* due to interferers lying in A n is represented by 



Ye'J>nA n 



P Y G T * Y \Y-T*\ 



(15) 



Therefore, I cc in (14) can be decomposed as I cc = Yln°=i ^i- 
To facilitate analysis, define a compound Poisson rv Z n as 



Z n — ^ PyGu'l 



(16) 



Ye<s>nA„ 



where {P Y Gu* Y } are i.i.d. and the number of terms in the 
summation, namely |$ n A n \, is a Poisson rv with mean I. 
Note that the distribution of Z n is independent of U*. Based 
on the geometry of A n , it can be obtained from (15) that 
In < {Vnp)~ a Z n . Since I = 2y/3p 2 \, it follows that 



(17) 



Combining (13) and (17) yields a lower bound on <p cc : 



^-^(g-^ip^'J- <18) 

Next, an asymptotic lower bound on f cc (l) as I — > oo can 
be derived by analyzing the large deviation of the summation 



in (18) as follows. As Z n is a sum over the i.i.d. sequence 
{PyGt'Y I Y G $ n A n }, it is necessary to characterize the 
large deviation of PyGjj-y as follows. 

Lemma 2. For sparse scattering and an arbitrary BS Y e 
$nC(r,p), E[PyGu, y ] is finite and 



-logPr(P F G(7.y > x) — 7rA 



7W 



x — > oo. 



(19) 



The proof of Lemma 2 is given in Appendix A. Analyzing 
the large deviation of Z n also requires the following result 
from [34, Proposition 7.1]. 

Lemma 3 ([34]). Consider a compound Poisson rv Z = 
Sm=i H"i where F follows the Poisson distribution and 
{H rn } are i.i.d. rvs independent of F. If the distribution of 
H m is either RV(r) with r > or WE(r) with < r < 0.5, 

Pi(Z - E[Z ] > x) ~ E[F] Pr( J ffi > x), E[F] -> oo 

if x > aE[Z a ] for all a > 0, where E[Z ] = E[F]E[i?i]. 

Since {PyG^y} G WE(t) with < r < 0.5 from 
Lemma 2, using the definition of Z n in (16) and applying 
Lemma 3 lead to the following result that is proved in 
Appendix B. 
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Lemma 4. Given x > 0, if a > 4 

-logPr(Z„ >fix) -ttA 
and if 2 < a < 4, 

-logPr(Z„ > ^i) y ttA I — 

V7^ 



e^Vx~, 



oo, 
(20) 

> oo. 
(21) 
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where n = 1, 2, • • • . 

Given Lemma 4, the application of the contraction principle 
from large-deviation theory (see e.g., [31, Theorem 4.2.1]) 
yields 4 



logPr (^Tn-'tZn > & 



-\ogVr(Z n >£^x) (22) 



\n=l 



as £ — > oo. Combining (18), (22) and Lemma 4 leads to an 
asymptotic lower bound on ip cc as shown below. 

Lemma 5. As I — > oo, the OPE for a cluster-center mobile 
satisfies 



ci^, a > 4 



c 2 ^t, 2<a<4 
where the constants c\ and c 2 are defined as 



(23) 



Cl 



2-y/3 \6l 



w^V^(2v / 3)^7" 



2) Asymptotic Upper Bound on the OPE: The OPE </? cc 
can be upper bounded by considering only the interferers for 
T* from a subset of the interference zone C(T*,p). For this 
purpose, define a "narrow" hexagonal ring 



A e =C(T*,y/T+ep)\C(T*,p) 



with e > and 



Z e — PyGt»y- 
re*n.4 e 



(24) 



(25) 



Note that Z f is a compound Poisson rv where the Poisson 
distribution has mean e£. Since \Y — T*\ < (1 + e)p for all 
y e $n A and A e e C{T*,p), I cc in (14) is lower bounded 
as 



I cc > [(l + e)p]- a Z e . 



(26) 



By combining (13) and (26) and using p = ^p, the OPE ip cc 
can be upper bounded as 

^ cc W<-logPr^-|(l + e)p) Z e >6- l wj. (27) 

Analyzing the scaling of the right-hand side of (27) as I — > 
oo leads to an asymptotic upper bound on ip cc as shown in 
Lemma 6 that is proved in Appendix C. 

Lemma 6. For sparse scattering and as £ — > oo, the OPE for 
a cluster-center mobile satisfies 

^ ^ (28) 
where the constant c\ is as defined in Lemma 5. 

4 The procedure is similar to that for obtaining (54) in Appendix B. 



3) Main Result and Remarks: Combining Lemma 5 and 6 
leads to the following theorem. 



Theorem 1. For sparse scattering and as I 
for a cluster-center mobile satisfies 

1) for a > 4, 



4ci 



oo, the OPE 



(29) 



2) and for 2 < a < 4, 



c 2 £f ^ ^ cc (^) ^ ^l, (30) 



where Ci and c 2 are as defined in Lemma 5. 
Several remarks are in order. 

1) Theorem 1 shows that <p cc (l) scales linearly with in- 
creasing i for a large path-loss exponent (a > 4) and 
at least sub-linearly for a moderate-to-small exponent 
(2 < a < 4). These results suggest that as I — > oo, P out 
diminishes exponentially and at least sub-exponentially 
for a > 4 and 2 < a < 4, respectively. The scal- 
ing of P out depends on a because it determines the 
level of spatial separation. Note that for a > 4, the 
asymptotic bounds on ip cc (£) have a ratio of 4/3 [see 
(29)] and hence are tight. Mathematically, the tightness 
of the bounds is due to the product rv PyGt*y in 
the expression for I cc in (14) having a distribution 
with a sufficiently heavy right tail, allowing accurate 
characterization of the asymptotic tail probability of I cc . 
However, as a decreases, the tail probability of PyGr-y 
reduces. This results in that the ratio of the asymptotic 
bounds on Lp cc (£) in (30) diverges as £ increases. 

2) It can be observed from Theorem 1 and the definitions 
of ci and c 2 that larger <p cc (£) results from increasing 
the ratio 5/-f, namely the minimum ratio between the 
magnitudes of beam main-lobes and side-lobes. In other 
words, as £ — > oo, the outage probability diminishes 
faster for sharper beams, agreeing with intuition. 

3) Theorem 1 suggests that for fixed outage probability, the 
outage threshold 9 should be proportional to 0. Cor- 
respondingly, the throughput of a cluster-center mobile, 
defined as R = log(l + 6), grows with increasing £ as 



R 



log t 



00. 



Note that R scales linearly with a because large a 
corresponds to more severe attenuation of inter-cluster 
interference. 



B. OPE for Typical Mobiles 

Consider the typical mobile U* and the corresponding OPE 
ip as given in (3). The asymptotic bounds on ip are derived in 
the following subsections. 
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1) Asymptotic Lower Bound on the OPE: First, a lower 
bound on ip is obtained as follows. For ease of notation, the 
distance Ly between the typical mobile and BS is re-denoted 
as L. As illustrated in Fig. 2(b), the interfering BSs for U* 
are Poisson distributed in the region A = C(T* , p)CMD{U* , L) 
where 0(A,r) represents a disk centered at A £ R 2 and 
with the radius r > 0, namely that 0(A, r) = {X e R 2 | 
\X - A\ < r}. Note that 0(U*,L) encloses the complete 
interference zone for U* due to the fact that any interfering 
BS for U* is farther than the serving BS Y* at a distance of 
L from [/*. As a result, the interference power for U* can be 
written as 

1 = J2 p yGu*y\Y -U*\- a 
Ye$nA 

< Yl p yGu>y\Y -U*\- a + 
Ye*nA 

]T P Y Gu* Y [m3x(\Y -U*\,L)]- a (31) 

FG#nC(T*,p)nO(T*,p) 

Using the facts that C(T* , p) = AU [C{T* ,p)nO(T*, p)] and 
\Y-U*\= max(|F - U*\, L) if Y G $ H A [see Fig. 2(b)], 
it follows from (31): 

/< V P F G a «y [max(|r-[/*|,L)]- Q (32) 
rG$nC(T*,p) 

oo 

<J2J2 PyGu-y [mnx(\Y-Y*\-L,L)]- a (33) 

n=l YGA n 

(33) is obtained based on the triangular inequality \Y— U*\ > 
|F - Y*\ - \Y* - and C(T*,p) = U„A with {A n } 
being the hexagonal rings defined in Section III-A1. Let D 
denote the distance from Y* to the boundary of C(T*,p): 



D = min Xe( j( T , p ) \Y* — X\, By the stationarity of the mobile 
and BS processes, Y* is uniformly distributed in C(T*,p), 
resulting in the following distribution of D: 

Pr(£> < x) = 1 - (l - , < x < p. (34) 

Since the shortest distance between Y* and a point in A n is 
\fnp — p + D, it follows from (33) that 

oo 

I< Y [™*{Vnp- P + D - L,L)]~ a Z n (35) 

n=l 

where Z n is defined in (16). From (3) and (35), 99 can be 
lower bounded as 

(00 
max( v / np - p+ 

n=1 v (36) 

D - L,L))~ a Z n > e^oj). 

Next, an asymptotic lower bound on 1^2 (£) is derived by 
analyzing the scaling of the right-hand size of (36) as I — > 00. 
For this purpose, it is shown in the following lemma that 
ip(£) can be asymptotically upper bounded by an expression 
comprising a series of the i.i.d. compound Poisson rvs {Z n }, 
which facilitates a similar approach as used for obtaining 
Lemma 5. 

Lemma 7. For sparse scattering and as £ —> 00, the OPE for 
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a typical mobile satisfies 



tp(£) y min maxmin 
^ w - \ z>0 



logPr [J2 n " 2Z n> 



\n=l 



logPr \ D a < y 



logPr ( L > — 




The proof of Lemma 7 is provided in Appendix D. By 
analyzing the scalings of the three terms in the lower bound 
on ip, an asymptotic lower bound on the OPE is obtained as 
follows. 

Lemma 8. For sparse scattering and as I — > oo, the OPE for 
a typical mobile satisfies 



1 - - ) log I 

a 



(37) 



The proof of Lemma 8 is provided in Appendix E. 

2) Asymptotic Upper Bound on the OPE: The analytical 
technique for deriving an upper bound on p essentially con- 
siders only interference to U* from interferers lying in a subset 
of the interference zone A defined in the preceding section. 
Specifically, define a region McA (see Fig. 3) as 



H 



= [A e C(Y*,D + L) n C(Y*, y/(D + L) 2 + 



--<Z(X-Y*)-Z(J Y . 



u) 
Y*) < 



-} 
6 J 



where u > and J Y * is a point in C(T* , p) such that | Jy* — 
Y*\ = D. Note that the hexagons in the definition of H are 
chosen such that the area of "H is a constant u/y/3. Then the 
OPE ip in (3) can be upper bounded as 



p(£) < - log (Pr (j2 Yen PyGu*y\Y -U*\- a > 

_1 w |$nw^uj Pr($ n H ^ 0) 



(38) 



Let Y denote an arbitrary BS in H conditioned on QDW ^ 0. 
Since |F -U*\ < \Y - Y* \ + L by the triangular inequality 
and 



v/3 



\Y -Y*\ < ^^(D + Lf+u 



(39) 



from the geometry of H (see Fig. 3), it follows from (38) that 



<p{l) <- logPr ( P Yo Gu» Yo > - Q [^^(D + LY + 



9 \ 2 



$ n H ^ - logPr($ n H i= 0). 



(40) 



By inspecting the scalings of the two terms at the right-hand of 
(40) as t — > oo, an asymptotic upper bound on ip(£) is obtained 
as shown in Lemma 9, which is proved in Appendix F. 

Lemma 9. For sparse scattering and as t — > oo, the OPE for 
a typical mobile satisfies 



fit) ■< \^gt 



(41) 



3) Main Result and Remarks: The following theorem re- 
sults from combining Lemma 8 and Lemma 9. 

Theorem 2. For sparse scattering and as i — > oo, the OPE ip 
for a typical mobile satisfies 



a I 2, 



(42) 



Several remarks can be made. 

1) The scaling of the OPE ip in Theorem 2 is largely 
determined by the left-tail probability [see (36) and (40)] 
of the distance D from the typical BS to the boundary of 
the affiliated cluster. The dominance of D in determining 
ip is due to that its distribution has a linear left tail [see 
(34)] that is heavier than the distribution tails of other 
random network parameters. As can be observed from 
(42), the asymptotic bounds on p are tighter for larger a. 
The reason is that the right tail of the interference-power 
distribution becomes lighter (with steeper slope) as a 
increases, which strengthens the mentioned dominance 
of D and thereby tightens bounds on tp. 

2) Theorem 2 shows that <p(£) scales logarithmically with 
increasing I. In contrast, from Theorem 1, the scaling of 
p cc (l) for a cluster-center mobile is much faster, namely 
at least sub-linearly with increasing I. The reason for this 
difference in the OPE scaling is that the typical mobile 
accounts for not only cluster-interior mobiles but also 
cluster-edge mobiles that are exposed to strong interfer- 
ence and as a result has much higher outage probability 
than a cluster-interior mobile. This suggests that cluster- 
edge mobiles are the bottleneck of network coverage and 
should be protected from strong inter-cluster interference 
by e.g., applying fractional frequency reuse [20] along 
cluster edges. 

3) The OPE scaling in Theorem 2 is closely related to 
the fact that the fraction of mobiles that are near 
cluster edges is approximately proportional to p^ 1 or 
equivalently £~2. Given the dominance of the outage 
probabilities for the cluster-edge mobiles over those 
of the cluster-interior mobiles, the outage probability 
for the typical mobile is expected to be approximately 
proportional to the fraction of cluster-edge mobiles and 
hence £~5. Consequently, the resultant OPE should be 
proportional to \ log I, which matches the result in 
Theorem 2. 

4) Unlike Theorem 1 (see Remark 3), Theorem 2 does not 
reveal the throughput scaling for a typical mobile. The 
reason is that the distribution of the distance D from 
a typical mobile to the boundary of the corresponding 
cluster [see (34)] dominates the OPE but is independent 
with the outage threshold 9 that determines the through- 
put. 

IV. OPE with Rich Scattering 

Sparse scattering is assumed in the preceding section. In 
this section, rich scattering is considered and the correspond- 
ing OPE is analyzed for cluster-center and typical mobiles 
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separately. It is shown that rich scattering decreases the OPE 
for cluster-center mobiles but has no effect on the OPE for 
the typical mobiles. 

A. OPE for Cluster-Center Mobiles 

1 ) Asymptotic Lower Bound on the OPE: The presence of 
rich scattering results in channel fading and hence affects the 
OPE. In particular, the resultant distributions of transmission 
power given channel inversion and interference-channel gains 
are characterized in Lemma 10 in the sequel. The effect of 
rich scattering is reflected in the difference between Lemma 2 
and Lemma 10. 

Lemma 10. For rich scattering and an arbitrary BS Y € $ f~l 

C(T*,p), as x — > oo, 

P.M,>.)~"3^.- (43) 

where T(-) denotes the Gamma function. 

The proof of Lemma 10 is given in Appendix G. Consider 
the lower bound on the OPE in (18) based on the sequence of 
compound Poisson rvs {Z n }, which also holds for ip cc with 
rich scattering. To analyze the scaling of the lower bound as 
t — > oo, the large deviation of Z n is characterized as follows. 

Lemma 11. For rich scattering and as I — > oo, 

> ^ _Lg-^+i (44) 



Pi(Z n >£^x) 
with n = 1, 2, ■ ■ ■ . 



(tt\)^x u 



The proof of Lemma 1 1 can be straightforwardly modified 
from that of Lemma 4 by applying Lemma 10 in place of 
Lemma 2; the details are omitted for brevity. It can be observed 
from (44) that the distribution of Z n does not have a sub- 
exponential tail as for the case with sparse scattering. This 
makes it difficult to apply the contraction principle as before to 
derive the scaling of the series X^^Li n ~^ Z n , which is needed 
for obtaining an asymptotic lower bound on ip. To overcome 
this difficulty, the current analysis applies the following result 
from [40, Theorem 2.3]. 

Lemma 12 ([40]). Consider a sequence of i.i.d. rvs {Z n } 
whose distribution belongs to RV(r) with r > and a 
sequence of nonnegative scalars {p n } with Y^n°=i Pn being 
finite for some < v < min(l,r). The tail probability of 
Y,n=i PnZn scales as 

/ OO \ OO 

Pr P^n > X ~ Pn Pr (^« >x), X^ OO. (45) 



\n=l 



n=l 



Based on Lemma 11 and Lemma 12, it is proved in 
Appendix H that as £ — >• oo, the OPE can be upper bounded 
as shown in the following lemma. 



Lemma 13. For rich scattering and as 
a cluster-center mobile satisfies 

1 



oo, the OPE for 



^> cc (£)h ^-l log*. 



2) Asymptotic Upper Bound on the OPE: The following 
lemma is proved using Lemma 10 and applying a procedure 
similar to that for proving Lemma 6 with the details omitted 
to keep the exposition precise. 

Lemma 14. For rich scattering and as £ — > oo, the OPE for 
a cluster-center mobile satisfies 



(47) 



3) Main Result and Remarks: The following theorem fol- 
lows directly from Lemma 13 and 14. 



Theorem 3. For rich scattering and as 
a cluster-center mobile satisfies 



oo, the OPE for 



^av - 1 j log* r< <p cc (l) r< ^av\og£ 



where v is the minimum signal diversity order. 



(48) 



(46) 



A few remarks are in order. 

1) By comparing Theorem 3 with Theorem 1, one can see 
that channel fading caused by rich scattering degrades 
ip cc dramatically. To be specific, as £ — >■ oo, ip cc (£) can 
scale at least sub-linearly with £ for sparse scattering 
but only logarithmically for rich scattering. Roughly 
speaking, fading increases the randomness in interfer- 
ence and thereby reduces the level of spatial separation. 
This introduces a larger number of significant interferers 
for the cluster-center mobiles with respect to the case of 
no fading and hence compromises the effectiveness of 
MCC. This is the key reason for the slower OPE scaling 
in Theorem 3 compared with that in Theorem 1. 

2) For single-cell transmissions over fading channels, in- 
creasing the BS density does not change the outage 
probability for an interference-limited network, as shown 
in [12]. In contrast, Theorem 3 indicates that it is 
possible to reduce outage probability by deploying more 
BS so long as the numbers of cooperating BSs increase 
proportionally. 

3) It is well-known that the effect of fading can be alle- 
viated by diversity techniques [41]. This is reflected in 
Theorem 3 where ip cc is observed to increase approxi- 
mately linearly with the minimum diversity order v if 
av is large. For this case, the asymptotic bounds on 
ip cc are observed to be tight. Moreover, ip cc also grows 
approximately proportionally with increasing a as inter- 
cluster interference is more severely attenuated. 

B. OPE for Typical Mobiles 

The type of scattering has no effect on the scaling of OPE 
for a typical mobile as £ — > oo as stated in the following 
theorem. 

Theorem 4. For rich scattering and as £ — > oo, the OPE for 
the typical mobile scales as shown in Theorem 2. 

The proof of Theorem 4 can be easily modified from 
that of Theorem 2 based on the new distribution of the rvs 
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{PyGwy} in Lemma 10. The detailed proof of Theorem 4 
is omitted. 

The insensitivity of ip with respect to the change on the 
scattering environment is due to that the distribution of D is 
independent with scattering and has a dominant effect on ip 
compared with the distributions of other network parameters 
(see Remark 1 on Theorem 2). Furthermore, since the distri- 
bution function of D is also independent with the diversity 
order N, it can be observed by comparing Theorem 3 and 4 
that unlike a cluster-center mobile, a typical mobile does not 
benefit from transmit diversity for improving the OPE scaling. 
Therefore, the result in Theorem 4 reiterates the importance of 
suppressing inter-cluster interference for cluster-edge mobiles 
to improve network coverage via MCC. 

V. Simulation Results 

The simulation method and settings are summarized as 
follows. The infinite network region is approximated by a disk 
centered at the origin, where BSs are Poisson distributed with 
density A = 10~ 2 and the disk area is chosen such that the 
expected number of BSs in the disk is 200 i.e., the disk area 
is 200/ A = 2 x 10 4 . The typical cluster region is centered 
at the origin and the size is determined by the expected BS- 
cluster size £. The main and side lobes of beams are uniformly 
distributed in the intervals [5,5'] = [6,10] and [0,7] = [0,1], 
respectively. Other parameters are sets as a — 4, N — 3 (for 
rich scattering), and 8 = 3. 

In Fig. 4, outage probability is plotted against increasing 
£ for different combinations of sparse/rich scattering and 
a cluster-center/typical mobile. To evaluate the asymptotic 
results derived in the preceding sections, Fig. 4 also displays 
curves obtained from the asymptotic bounds on the OPE 
as follows. Consider a typical mobile and let cp + and 
represent the asymptotic upper and lower bounds on the OPE, 
respectively. Note that outage probability can be approximated 
as P ut ~ be~ v ^ if £ ^> 1 where b is a constant. For this 
reason, the functions b\ e~ v W and b 2 e~' p ^ are plotted in 
Fig. 4 and identified by the legends "OPE asymptotic upper 
bound" and "OPE asymptotic lower bound", respectively, 
where the constants bi and b 2 are chosen such that the 
matching analytical and simulation curves overlap at their 
rightmost points for ease of comparison. Similar curves are 
also plotted in Fig. 4 for a cluster-center mobile. The curves 
based on analysis and simulation are observed to be closely 
aligned if £ is sufficiently large, indicating that the derived 
asymptotic bounds on the OPE (especially the asymptotic 
lower bound) are accurate. In particular, for the cluster-center 
mobile with sparse scattering, the curve from the asymptotic 
lower bound on the OPE overlaps with the simulation curve 
and hence this bound is tight even for small values of £. 

Next, it can be observed from Fig. 4 that as £ increases, 
the outage probability for a cluster-center mobile decreases 
rapidly but the outage probability for a typical mobile remains 
almost unchanged and close to the result for the case of 
no MCC (specified in Fig. 4 using dashed lines). In other 
words, it is verified that MCC benefits only cluster-interior 
mobiles and cluster-edge mobiles limit network coverage. This 



observation is consistent with findings from implementing 
MCC in practical networks [19], [42], [43]. Furthermore, with 
respect to sparse scattering, rich scattering is observed to 
increase outage probability for cluster-center mobiles by up 
to several orders of magnitude. 

Fig. 5 compares the outage capacity of cluster-center and 
typical mobiles, namely their maximum throughput given the 
maximum outage probability of 0.05. The observations from 
Fig. 5 agree with those from Fig. 4. Specifically, the outage 
capacity for a typical mobile is marginal even as £ increases 
while the capacity for mobiles without MCC is approximately 
zero. In contrast, the outage capacity for a cluster-center 
mobile increases rapidly with growing £ and sparse scattering 
results in much higher capacity than rich scattering. 

VI. Concluding Remarks 

In this paper, a novel model was proposed for a cellular 
downlink network with MCC. The network coverage was 
analyzed in terms of the outage-probability exponent. It was 
shown that though the performance gain for cluster-interior 
mobiles from MCC is large, the gain for a typical mobile is 
small as it is likely to be located near the edge of a base-station 
cluster and exposed to strong inter-cluster interference. This 
finding provides an explanation for the marginal gain of MCC 
in practice, and suggests the need to design a new medium- 
access-control protocol or apply fractional-frequency reuse for 
protecting cluster-edge mobiles. 

This work opens several interesting directions for further 
research. In particular, instead of using a lattice, base-stations 
can be clustered by a random process such as a Poisson 
random tessellation that gives non-uniform expected BS- 
cluster sizes. Moreover, the current interference-coordination 
algorithm that requires multi-antennas at base stations can 
be replaced with a network-MIMO algorithm that supports 
cooperation between single-antenna base stations at the cost 
of inter-cell data exchange. Last, the proposed analytical 
framework can be applied to study the performance of other 
MCC algorithms and heterogeneous networks with MCC. 

Appendix A 
Proof of Lemma 2 

Consider an arbitrary BS Y £ For convenience, define 
j3 = Wy /Gjj*y with support [S/j,oo) and the probability 
density function is denoted as fp. Using Py = ujL y /Wy 
from (12), 

Pr(P r GW >x)= / Pr (L Y > ™) fp{r)dr, x > 0. 

7 

Substituting the distribution function of Ly in (5) gives 

e -*M™/«)« f ( T )dT. (49) 

7 

The right-hand size of (49) can be expanded for e > as 



Pt(P y G 



u*y >x) = 



t\{tx/u)) t 



U(r)dr+ 



s: 



e -**(™/u>)<« f p ( T )dT. 



(50) 
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(a) Sparse scattering 
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(b) Rich scattering 

Fig. 4. Outage probability versus expected BS-cluster size i with (a) sparse scattering or (b) rich scattering. For comparison, the outage probability for the 
case of no MCC is specified by dashed horizontal lines. 



It follows that Next, it can be obtained form (50) that 



Pr(PyG t /.y > X 



) > e -^(( 1 + £ )^)" x Pr(P y GW > x) < e- A (( 1+£ )^)° Pr U > (1 + e)^\ + 



'Y '-yy' e --^^r p r _ < < (1 + e ) 

7 V7 7 

Thus, as x -> oo, As a result, 

-logPr(PyGr/*y > x) < (1 + e)-vrA ( — V x« . (51) - logPr(PyGr/»y > x) h ttA ( — ) x« , a: -> oo. (52) 
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Fig. 5. Outage capacity per mobile versus expected BS-cluster size I for the maximum outage probability of 0.05. Note that the outage capacity per mobile 
for the case of no MCC is approximately zero. 



Combining (52) and (51) and letting e — > yield the desired 
result in (19). Last, the claim of E[PyGjj*y] being finite 
follows from the fact that E[L a ], E[Gu* Y ] and EfW^ 1 ] are 
all bounded, which follows from the distribution of L in (5) 
and those of G\j*y and Wy in Lemma 1. ■ 

Appendix B 
Proof of Lemma 4 

First, consider the case of a > 4. Since PyGjj*y € WE(r) 
with < r < 0.5 according to Lemma 2 and E[Z n ] = 
£E[PyGu*y]> applying Lemma 3 gives the desired result in 
(20). 

Next, consider the case of 2 < a < 4. It is claimed that as 

I — !• oo, 

-logPr(Z„ >£%x)h 

' /r , „ , 4(1+g) _ „«_\ (53) 



X — > oo, 



logPrl {PyGu^y)^ >£%x 



with e.x > 0. To prove this claim, let V\, V%, ■ ■ ■ , Vf. denote 
k i.i.d. rvs following the same distribution as PyGjj*Y for an 
arbitrary Y G $ (1 C(T*,p), By using Lemma 2 and applying 
the contraction principle from large-deviation theory [31, The- 
orem 4.2.1], for a set of nonnegative numbers {xi,X2, ■ ■ ■ , xt} 
and as x — >• oo, 



logPr iJ2 V k > x 

\n=l j 



ttA 



6 

7W 



inf 



(54) 



where (54) results from the inequality ^Xm=i — 
Sn=i a?? if < p < 1. It follows from Lemma 2 that as 



•logPr(Vi " >xj-7rAf — ) - - . (55) 



Using (55) and again applying the contraction principle give 



logPr [Y,V k " >x UttA 



70; 



a;2(i+«) (56) 



as x — y oo. Given 2 < a < 4, comparing (54) and (56) yields 
- log Pr (j2 V >> > ^ - log Pr ^ > (57) 



as x — ?> oo. Since the inequality in (57) holds for arbitrary k, 
the claimed inequality in (53) is proved. Recall that PyGwY 
for arbitrary Y 6 $ fl C(T*,p) has the same distribution as 
Vi. By inspecting (55), P Y Gu* Y € WE(r) with < r < 0.5. 
Therefore, it can be derived similarly as (20) in the lemma 
statement that 



logPrf ]T (P Y Gu* 
\Ye<s>nA n 



, *(i+<0 „a 
y) « >£ 2 X 



(58) 



7rA I — ) f'(i+«)j;2P+<) . 



Substituting (58) into (53) and letting e — > gives (21) in the 
lemma statement. ■ 
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Appendix C 
Proof of Lemma 6 

From the definition of Z t in (25) and (27), 
^(£) < - log (Pr + e )^ P r GW > - Q 



(59) 



Since Pr($ ni f ^) = (l- e - £f ) and I = 2^/Zp 2 \, using 
(59) and applying Lemma 2 give that 

The desired result follows from (60) by letting e — > 0. ■ 

Appendix D 
Proof of Lemma 7 

By expanding outage probability defined in (1), 
Pout <Pr(l>^\L< I pJ+Pr(L> 1 pJ. (61) 



The substitution of (35) yields that 



Pout < Pr ( J2 Zn ( max (V™P -p + D-L,L)) a >"r 



, n=l 



X < 



Pr L > 



< Pr ( ^ Z n ( V^p - p + ^ ] > " | + 



\n=l 



(62) 



Since D < p, the replacement of p in (62) with D/2 further 
upper bounds P out as 



Pout<Pr(^n-^„> — 



ri /• " ) (63) 



Applying the similar method as for obtaining (62) results in an 
upper bound on the first term on the right-hand side of (63): 



Pr ( " 1 inin 



\n=l 



2 Q 6> 



z>0 



Pr J^n-f Z„ > 

\n=l 



2 a ( 



+ Pr(D a < z) 



By combining (63) and the last inequality, 



Pout < min 

z>0 



V7l=l 



LOZ 



+ Pr(L> Q < z) 



Pr,L> ? 



(64) 



As f — » oo, the term on the right-hand size of (64) that 
decays at the slowest rate dominates the other two terms. 
Specifically, given (64) and the definition of ip in (2), applying 
[31, Lemma 1.2.15] yields the desired result in the lemma 
statement. ■ 



Appendix E 
Proof of Lemma 8 

Consider the three terms in the asymptotic lower bound on 
ip in Lemma 7. By setting 2~ a 0- 1 ujz = £[Z n ]+ p 2{1+ ^ with 
e > 0, the procedure similar to that for obtaining Lemma 4 
can be applied to derive the following asymptotic lower bound 
on the first term: as I — > oo, 



-logPr(5>-fZ„>H)b 



, n=l 



(65) 



ttA 



7W 



2\/3A 



The scaling of the second term is obtained using (34) and the 
aforementioned constraint on y as 



- logPrp" < z )~Ul- 2 i 1 ±^ ) loo , 



(66) 



as I — > oo. Using the distributions of L and D in (5) and (34) 
respectively, 



Pr [L > 



D\ 2 f p 

2) = - P L ' 



1 - - dT 
P) 

p — » oo. 



VXp ^\p 2 ' 
By substituting I = 2\/3p 2 , the third term scales as 



- log Pr L > 



D 



■tog*, 



oo. 



(67) 



Last, the substitution of (65), (66) and (67) into the asymptotic 
lower bound on <p> in Lemma 7 and letting e — > lead to the 
result in the lemma statement. ■ 

Appendix F 
Proof of Lemma 9 

As the area of % is u/\/3, the number of BSs in %, namly 
|3> fl H\, follows the Poisson distribution with mean uX/y/3. 
Using this fact and (40), 

<fW < - log ^Pr {py G v . Yo > | (^^(D + L) 2 +u+ 
L ] |£> < a;] Pr(£> < a;)) -log(l-e"^) 



<-togPr fV GW„ > 



y / {x + L) 2 + U+ 



Lj J -logPr(L> < x) -log(l-e vs) 

with x > 0. By keeping u and x constant and letting p — > oo, 

¥>W ^ -logPr(D < x) 

~logp (68) 

where (68) uses the distribution of Z) in (34). The substitution 
of I = 2\fi\p 2 into (68) proves the desired result. ■ 
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Appendix G 
Proof of Lemma 10 

Consider an arbitrary BS Y e $ and the corresponding 
parameters {Py,Gy,Wy,Ly}- The subscripts of these pa- 
rameters are omitted in the remainder of the proof to simplify 
notation. Given P = ujL a /W from channel inversion and 
x > 0, it follows from Lemma 1 that 



Pr(PG > x) = Pr (G > u~ 1 WL" 'x) 

-oj- 1 WL- a x~ 



= E 

= E 

= E 

> E 



Jo 



N-l 



(N-iy. 



e~ T dr 



1 



(1 + u-T-xL-t) 1 * 

1 X 

(1 +u- 1 xL~ a ) N ' uL c 

(X 



(69) 
(70) 
(71) 



> log a: 



Pr 



\ujL c 



> logx^ 

u N L aN x 
(1 _|_ ^jfL)N x N I ~ZAJ* 

Pr ( — %— > logx) 



> logx 



> E 



N t qJV 



uj n L 



> logx 



= E 



Pr (^ >l0gX ) 



N 



,,N 



(72) 



r(^ + i) 



-N 



.N 



N T(<f + 1) Pr(/V = i/) 
(7rA)^ 



X — > oo 



(73) 



where (69), (70) and (72) use the distributions of G and W 
in Lemma 1 and L in (5), respectively. Next, from (71), 



.N-N TCtNl 



Pv(PG > x) < E [uj"x-"L° 



w N T (^f + 1) Pr(iV = i/) 
(ttA)t 



(74) 



as x — > oo, where (74) is obtained using a similar procedure 
as (73). Combining (73) and (74) gives the desired result. ■ 

Appendix H 
Proof of Lemma 13 

To apply Lemma 12, define Z' n — Z n /i^ and z — 0~^x 
and rewrite (44) as 

, .^r (f + 1) Pr(iV = v) _ v 



Pr(Z' n > z) 



(irX)'- 



(75) 



as z -> oo. It can be observed from (75) that Z' n e RV(^). 
Moreover, given a > 2, the sum Yl^Li n ^ can checked 



to be finite. Therefore, using (75) and applying Lemma 12, 



pr(5>-*z;> 



\n=l 



w ^r + 1) Pr(JV = i/) 

" / x 

(7rA) 2 

oo 

n=l 

Combining the definitions of Z' n and z and (76) yields 
Prf£n-fZ„-^-^ r ^ + 1 ) P ^ = ^ 



(76) 



> € 2 x 



(TTA)Tf 

OO 



(77) 



The desired asymptotic lower bound on the OPE follows from 
(18) and (77). ■ 
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